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We investigate the electronic state and structure transition of BaNi2As2, which shows a sim- 
ilar superconducting phase diagram as Fe-based superconductors. We construct the ten-orbital 
tight-binding model for BaNi2As2 by using the maximally localized Wannier function method. The 
Coulomb and quadrupole-quadrupole interactions are treated within the random-phase approxima- 
tion. We find that the charge susceptibility of quadrupole Ox^-y 2 nas the largest peak structure 
at q = (tt, 0, 7r) which is derived from both the in-plane and out-of-plane oscillations of Ni ions. We 
propose that the antiferro-quadrupole O x z-y 2 order with q = (jr, 0,7r) is the origin of the zigzag 
chain structure reported in experiments. The obtained antiferro-quadrupole fluctuations would then 
be the origin of the strong coupling superconductivity in Ni-based superconductors. We identify 
similarities and differences between Ni- and Fe-based superconductors. 

PACS numbers: 74.70. Xa, 75.25.Dk, 74.25. Kc 



I. INTRODUCTION 

Since the discovery of superconductivity in the Fe- 
based pnictide by Kamihara et a/., 1 great efforts have 
been devoted to uncovering its mechanism^ The un- 
doped Fe-based superconductors exhibit structural tran- 
sition from tetragonal to orthorhombic at T s and spin- 
density- wave (SDW) antiferromagnetic transitions at 
Tn^~— In many compounds, T s is higher than Tn, and 
T s is expected to be realized by orbital polarization^ - — 
In both electron- and hole-doped cases, superconduc- 
tivity appears near the structure, orbital, and SDW 
orders^"— i 12 i 13 It is, therefore, natural to consider the 
possibility that the orbital and spin fluctuations pro- 
vide the pairing mechanism, and two types of supercon- 
ducting states are proposed: s± superconductivity medi- 
ated by the spin fluctuation^ - — and s ++ superconduc- 
tivity mediated by the orbital fluctuation*^ - — In many 
compounds, the superconductivity is very robust against 
impurities j22 - — indicating that the s ++ -wave state is re- 
alized in these compounds.— 

BaNi2As2 also exhibits both structure and supercon- 
ducting transitions^ - — whereas no evidence of SDW 
order has been observed. Above the structure transi- 
tion temperature T s ~ 130 K, it has the same tetragonal 
structure as BaFe2As2, which is a typical parent com- 
pound of Fe-based superconductors. At T s , BaNi2As2 
exhibits a first-order triclinic structure transition ^ 26 ! 27 
whereas BaFe2As2 exhibits a second-order orthorhom- 
bic transition. 5,6 Below T s , the Ni atoms form zigzag 
chains with shorter Ni-Ni distances (~ 2.8 A) and the 
chains are separated by significantly longer Ni-Ni dis- 
tances (~ 3.1 A) #22 The superconductivity emerges at 
T c ~ 0.7 K in this triclinic phase, which is considered to 
be of the conventional weak-coupling BCS type j 26 i 28- — 
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From a previous study of density functional theory 
(DFT) calculations for BaNi2As2, it was found that a 
strong electron-phonon interaction (A ep = 0.76) is in- 
duced by the low-energy (~50 K) Ni and As oscillat- 
ing modes #22 Kudo et al. recently reported the suppres- 
sion of the structural transition in the P-doped com- 
pound BaNi2(Asi_ ;E P :E )2 and the emergence of a new 
strong-coupling superconducting state (AC/jT c ~ 2) 
with T c > 3.3 K for x > 0.07. 30 

In Fe-based superconductors, a large softening of the 
elastic constant Cqq has been reported above T s £^~— 
Also, spontaneous in-plane anisotropy of the electronic 
states, so called the nematic state, was reported by 
the magnetic torque^ and resistivity measurements ] 35 ! 36 
These facts indicate the existence of large orbital fluctu- 
ations, and the relationship between the orbital fluctu- 
ations and superconductivity has attracted much atten- 
tion. However, the softening and the relation T s > Tn 
can not be explained within the random phase approx- 
imation (RPA) for multi-orbital Hubbard model^ Re- 
cently, we have improved the RPA by including the vertex 
corrections (VCs), and found that the strong ferro- and 
antiferro-orbital fluctuations are caused by the strong 
spin-orbital coupling described by the VC. 8,9 The ob- 
tained orbital fluctuations well explain both the struc- 
ture transition, including the large Cqq softening, and 
s ++ -wave superconducting state. 

For BaNi2(Asi_ ;E P ;E )2, a large phonon softening to- 
ward the structure transition has also been reported, 30 
although the structure transition is of first order. There- 
fore, strong orbital fluctuations are expected to exist in 
BaNi2(Asi_ ;E P :E )2 as in Fe-based superconductors. It 
is noteworthy that lii- x Pt x Te2 37 and CaC^^ also 
exhibit superconducting transitions and a lowering of 
phonon frequencies next to the structure transitions. 
Therefore, it is important to understand the relationship 
among structure transition, orbital fluctuations, and su- 
perconductivity. It is also significant to identify similari- 
ties and differences between Ni- and Fe-based supercon- 
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ductors. 

Herein, we investigate the electronic state and struc- 
ture transition of BaNi2 As2 based on RPA. The obtained 
critical value of the Coulomb interaction for the SDW 
state, U c , is very large, which is consistent with the ab- 
sence of the SDW state. However, it is found that the 
quadrupole interaction resulting from Ni-ion oscillations 
gives rise to the antiferro-quadrupole (AFQ) order, which 
presents a natural explanation for the zigzag chain struc- 
ture in the triclinic phased We also discuss the orbital- 
fluctuation-driven superconductivity in Ni-based super- 
conductors. 



II. POSSIBLE QUADRUPOLE ORDER 

Before performing numerical calculations, we discuss 
the possible quadrupole order that represents the exper- 
imental zigzag chain structure. Figure QJa) and DJb) 
present schematic pictures of the Ox 2 -y 2 order with 
momentum q = (0,7r,7r) that corresponds to the zigzag 
chain structure. Here, a, 6, and c are lattice constants, 
and there are two Ni ions, A and B, per unit cell. In 
this study, we set the X and Y axes along the a and 
b axes (so that the Ni-As direction is in the a-b plane) 
and the Z axis is then parallel to the c axis, as shown in 
Fig-DJa). The symbols + and — around the Ni ions rep- 
resent the sign of the charge distribution of the Ox 2 -y 2 
quadrupole. Ox 2 -y 2 is an Ising-type order parameter. 

As shown in Fig.[TJb), the values of Ox 2 -y 2 on Ni ions 
have the same sign on each zigzag chain (as shown by 
dotted lines). In this case, the bond length on the zigzag 
chain would become shorter owing to the electrostatic po- 
tential. However, Ni ions on different zigzag chains would 
feel a repulsive force and form a longer bond. Experimen- 
tally, longer bond (shorter bond) is about 3.1 A (2.8 A)i22 
Such a large difference between longer and shorter bonds 
might indicates the importance of electron-phonon inter- 
action. Since Ox 2 -y 2 <x ^d X z ~ n d Y z^ the quadrupole or- 
der Ox 2 -y 2 > corresponds to the orbital polarization 
n^ xz > n^ YZ ^ Thus, the AF Ox 2 -y 2 order is equivalent 
to the antiferro orbital order shown in Fig.QJc), in which 
the mainly occupied d-orbital wave functions are shown. 

In this paper, we investigate the development of 
quadrupole fluctuations due to small quadrupole interac- 
tion, which is induced by Ni-ion oscillations. We obtain 
the divergent development of Ox 2 -y 2 susceptibility at 
q = (7r,0,7r) (and q = (0,7r,7r)), which is consistent with 
the zigzag structure shown in Fig. l(a)-(c). 



III. FIRST PRINCIPLE STUDY AND MODEL 
HAMILTONIAN 

First, we perform a DFT calculation for BaNi2As2 
with the generalized gradient approximation by using 
the wien2k packaged In our calculation, we used ex- 
perimental lattice parameters^ Figure [2](a) and Efb) 




FIG. 1. (Color online) (a) Crystal structure of BaNi2As2 
with AF X 2 _ Y '2 order with momentum q — (0,7r,7r). The 
dotted zigzag lines correspond to the shorter bond. A and B 
represent the Ni-A and Ni-B ions in a primitive unit cell. The 
box with a solid line corresponds to a conventional unit cell, 
(b) Two-dimensional picture of the AFQ. The square with 
a solid line corresponds to a primitive unit cell. The arrow 
(0, 2a) represents the translation vector corresponding to q = 
(0,7r). (c) The antiferro orbital order, which is equivalent to 
the AFQ order in (b). 



show the band dispersions obtained from the DFT cal- 
culation (solid lines) on the tetragonal and the triclinic 
phases, respectively. The Brillouin zone is shown in 
Fig.[2jc). The band structure of BaNi2As2 in the tetrag- 
onal phase is qualitatively similar to that for BaFe2As2, 
although the Fermi level is shifted upward since Ni 2+ 
contains two more valence electrons than Fe 2+ . The ob- 
tained result is consistent with those of previous studies 
of DFT calculations ! 29 ! 41 " — As a result of this upward 
shift, BaNi2As2 exhibits very different electronic prop- 
erties from those of BaFe2As2- The Fermi surfaces in 
tetragonal and triclinic phases are shown in Fig. Efd) and 
EJe), respectively. As shown in Fig. EJd), there are three 
large Fermi surfaces FS1-FS3, a small dishlike electron 
pocket, FS e , around Z point (0, 0, 27r), and a hole pocket 
FSh around (7r, 0, tt). FSi and FS2 have simple cylindrical 
forms, whereas FS3 has a very complex structure. The 
shape of these Fermi surfaces is similar in both tetrago- 
nal and triclinic phases. In contrast, FS e disappears in 
the triclinic phase, as shown in Fig. Efe). Further, FSh 
disappears around (±7r, 0, ±tt) owing to the energy gap 
opening up as the symmetry is lowered, and it becomes 
small at around (0, ±7r, ±7r). We note that the X and Y 
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FIG. 2. (Color online) (a) and (b) Band structure of 
BaNi2As2- In (a), the solid and dashed lines correspond to 
the DFT calculation for the tetragonal phase and the tight- 
binding model for AE = 0, respectively. In (b), the solid lines 
correspond to the DFT calculation for the triclinic phase, 
and dashed lines are given by the tight-binding model for 
AE — —0.5 eV. (c) Brillouin zone, (d) Fermi surfaces in the 
tetragonal phase given using DFT. (e) Fermi surfaces in the 
triclinic phase. 



axes are not equivalent in the triclinic phase. 

Next, we introduce the ten-orbital tight-binding model 
for BaNi2As2 with the tetragonal structure. The tight- 
binding Hamiltonian is given by 



>.,l,a 



i a,l,a 
i,j ct,/3 ,l,m,<j 




FIG. 3. (Color online) (a) The band structure of the tight- 
binding model for AE = and AE = -0.5 eV. (b) and (c) 
The Fermi surfaces of the tight-binding model for AE = and 
AE = — 0.5 eV projected onto the kx-ky plane. The arrows 
(7r,7r) and (7r,0) are the nesting vectors that correspond to 
the peak of the quadrupole susceptibilities shown in Fig. [5] 



where i and j denote the unit cell, a and f3 represent Ni 
ions A and B, I and m denote the d orbital (d 3 z 2 -R 2 , 
dxz, dyz, d X 2 -Y 2 i or dxy), and a = ±1 is the spin 
index. c^ x a (cL m a ) is a creation (annihilation) operator 
of the Ni 3d electron. The on-site energies e/ and hopping 
integrals t % £ a l m are obtained from maximally localized 
Wannier functions using the WANNIER90 code^ and the 
wien2wannier interfaced 

We make comparison of the band structures given 
by the DFT calculation and the present tight-binding 
model. In Fig. Efa), dashed lines represent the band 
dispersions of the obtained tight-binding model. Af- 
ter the triclinic structure transition, band dispersions 
show (i) energy level splitting with symmetry lowering. 
In addition, we found (ii) a large orbital-dependent en- 
ergy level shift without symmetry lowering, reflecting 
the first-order transition. We verified that change (i) is 
much smaller than change (ii). The latter change can be 
roughly reproduced by energy level lowering of dxz and 



dyz orbitals, AE = AEd 



AE d 



-0.5 eV. The 



obtained tight-binding dispersions for AE = —0.5 eV 
are shown in Fig. [2fb) by dashed lines. Note that the 
tetragonal symmetry is not violated by AE. 

For comparison, we show the band structure for AE = 
and AE = -0.5 eV in Fig. E^a). We also show 
the Fermi surfaces (projected onto the kx~ky plane) for 
AE = and AE = -0.5 eV in Fig.^b) and^c), respec- 
tively. In Fig.[3fc), FS e disappears and FSh shrinks, con- 
sistently with the DFT band dispersion for the triclinic 
phase. Hereafter, we calculate both cases AE = and 
AE = — 0.5 eV, and we discuss the spin and quadrupole 
fluctuations. 
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IV. NUMERICAL RESULTS BY RPA 
A. Formalism 

Here, we present the numerical results obtained by 
RPA. The noninteracting susceptibility for zero fre- 
quency is given by 



,m,m 



(a) = _1 v /M fc + <?))-/M fc )) 

xu a ,i,b(k + q)u* ^ b (k + q) 

XUp,m',b'(k)u* llb ,(k), 



(2) 

where e&(fc) and /(eb(fc)) are the dispersion and the Fermi 
distribution function of a quasiparticle for band b with 
momentum k. The function u a ^^(k) = (a,l\b,k) con- 
nects the orbital and the band spaces, given by the pro- 
cedure for diagonalization of H°. 

Now, we analyze the effects of the Coulomb and 
quadrupole-quadrupole interactions using the RPA. The 
on-site Coulomb interaction is composed of the intra- 
orbital term U, the inter-orbital one U' ', Hund's cou- 
pling J, and pair hopping J' . We set these parameters as 
U' = U — 2J and J' = J to satisfy the rotational invari- 
ance of the electron-electron interaction term. The bare 
four-point vertex for the spin channel is given by 
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(3) 



Further, the bare four-point vertex for the charge channel 
is given by 



-C — 2(y qua d) 



with 



U Z,Z', m,m' 



= < 



-U f + 2 J, 
2U f - J, 



I — V = m = m', 

I = m 7^ V = ra', 

I = V 7^ m = ra', 

I = m! = m, 



(4) 
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and 



(Vqua.d)i j' ?m , m ' 
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(6) 



where g and g' are the charge quadrupole-quadrupole 
coupling constants. In Eq. (j6j) the quadrupole inter- 
actions for channels XZ(Y Z) and X 2 — Y 2 are medi- 
ated by in-plane and out-of-plane oscillations of Ni ions, 
respectively The charge quadrupole operator is de- 

and the coeffi- 
) for T = XZ. 



>.,l,(j^a,m,(j'> 
I X Z 



ELm i f 

■xz^7(l\H\ m ) 
In the RPA framework, the spin (charge) susceptibility 

is obtained by using T s ^ as 46 



fined as Of ^ Lmf 
cient 6r is defined as o 



X s{c \q) = 



i_ r s(c)^o (g ) 



(7) 



(a) x sAA (q) 

AE=0 



(b) x sAA (q) 
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FIG. 4. (Color online) The obtained total spin susceptibility 
X sAA (q) for (a) AE = and (b) AE = -0.5 eV. 



The ordered state is realized when the Stoner fac- 
tor a s ( c ) = 1, which is the maximum eigenvalue of 
r s ( c ) x s ( c \q) . Hereafter, we set the band filling as n = 8.0 
and the temperature as T = 0.02 eV. We also set J = U/6 
and g = g r . 



B. Spin fluctuations 

First, we calculate the total spin susceptibility result- 
ing from to Coulomb interaction: 



X 



s a, (3 



(<z) 



s a,f3 

Xi 



(<?)• 



(8) 



Figure a) shows the obtained X sA ' A (#) f° r = 
and U = 2.32 eV; the corresponding spin Stoner fac- 
tor is a s = 0.98. X sA ' A ((?) nas its highest peak at the 
incommensurate wave vector q & (7r,0,0). Figure H^b) 
shows x sAA (q) for AE = -0.5 eV and U c = 2.07 eV 
(a s = 0.98). In both cases, X sA ' A ((?) nas a peak at q = 
(7r,0). It seems that the intra-orbital nestings of Fermi 
surfaces, which give the enhanced spin susceptibility, are 
not significantly different in the two cases. (In contrast, 
BaFe2As2 has spin fluctuations with q = (7r, 7r)^^ which 
corresponds to the stripe-type SDW order. 5 7 ) The ob- 
tained critical value U c > 2 eV is about twice as large 
as the value reported for iron-base superconductors us- 
ing RPA calculations j & 14 i 19 This fact implies that a large 
Coulomb interaction is required to drive the SDW tran- 
sition in BaNi2As2- We stress that the obtained U c is 
largely underestimated since the self-energy correction is 
neglected. Therefore, we conclude that the spin fluctua- 
tions are very small in BaNi2As2- 



C. Orbital fluctuations 

Next, we calculate the charge quadrupole susceptibili- 
ties resulting from the quadrupole interaction^ 



( \ m.m /r\\ 

A<l) r ■ ( 9 ) 
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FIG. 5. (Color online) The charge quadrupole susceptibility 
X? AA (q) for [(a) and (b)] AE = and [(c) and (d)] AE = 
—0.5 eV. The peak structure of Xx 2 -y 2 (q) f° r Ai£ — —0.5 eV 
in (d) is consistent with the experimental zigzag structure in 

Fig. in 



Figure [5fa) and Mb) show the obtained Xxz' (#) anc ^ 
Xx2 A -y 2 (g) for AE = 0, £7 = and g = 0.392 eV, re- 
spectively. The corresponding charge Stoner factor is 
a c = 0.98. When AE = 0, X X z{y Z )(<i) has its hi g n " 
est peak at q = (7r,7r,7r), which is given by the inter- 
orbital nesting between the d X z(YZ) orbital on FS2 and 
the dx 2 -Y 2 orbital on FS3. The nesting vector is shown 
as the thick arrow in Fig. [3fb). The obtained peak po- 
sition at q = (71", 7r, 7r), however, is inconsistent with the 
structural transition of BaNi2As2- 

We also show Xxz^iv) ano ^ Xx^-Y 2 (g) ^ or = 
-0.5 eV, U = and g = 0.399 eV in Fig. Efc) and 
Efd), respectively. The charge Stoner factor is a c = 0.98. 
Compared to the case of AE = 0, the peak of Xxz(yz) (#) 
at qr = (7r, 7r, 7t) is suppressed, and instead, Xx^-^y 2 (#) 
has its highest peak at qr = (7r,0,7r) (and </ = (0,7r,7r)) 
owing to the inter-orbital nesting between the d X 2 -Y 2 
orbital on FSi and the d 3Z 2 -R 2 orbital on FS2. The 
nesting vector is shown as the thick arrow in Fig. [3fc). 
The corresponding AF Ox 2 -y 2 order with q = (71", 0,7r) 
is consistent with the zigzag chain structure reported 
experimentally, 27 as already shown in Fig. HJ The ob- 
tained change in the inter-orbital nesting conditions orig- 
inates from the small modifications of FSi and FS2 by 
AE. Although FS e and FSh are mostly changed by AE, 
the contribution of these changes to the susceptibilities 
is negligible since these pockets do not contribute to the 
dominant nesting. 

The enhancement of y^ a ^(q) is mainly caused 



by the quadrupole interaction with respect to op in 

Eq. ©. Thus, XxzIyz)^) and Xx^-y 2 ^) are me_ 
diated by the in- and out-of-plane oscillations of Ni 
ions, respectively^^"— In more detail, both suscepti- 
bilities Xxz{yz)(^ ano - Xx^-y 2 ^) enhance each other 
at q = (7r, 0, 7r) through quadrupole off-diagonal term 
resulting from the orbital hybridization. 

V. DISCUSSIONS 

We have shown that experimental zigzag chain struc- 
ture is reproduced by the divergence of Xx^-y 2 ( 7r ' 0> 7r ) 
due to quadrupole interaction g. The obtained critical 
values g c ~ 0.4 eV for the quadrupole order are about 
twice as large as the values reported for Fe-based super- 
conductors (g^ e ~ 0.2 eV) using RPA calculations 
However, a strong electron-phonon interaction of A ep = 
0.76 for BaNi2As2 is found using a DFT calculation, 29 
whereas ~ 0.2 for LaFeAs0 . 15 i 47 Therefore, a charge 
quadrupole interaction mediated by an electron-phonon 
interaction would be dominant for BaNi2As2, whereas 
the Coulomb interaction would be less important. 

The result of the present analysis is consistent with 
the previous DFT study, which verified that the total 
energy is lowered by the formation of the zigzag chain 
structure.— In addition, the present study clarifies that 
the origin of the zigzag structure is the AF Ox 2 -y 2 order 
induced by an electron-phonon interaction and the inter- 
orbital nesting of the Fermi surfaces for AE = — 0.5 eV. 

Finally, we discuss the role of the AFQ fluctuations on 
the superconductivity. We have obtained two types of 
AFQ fluctuations, OxziYZ) with q = (tt, 7r, tt) for AE = 
and Ox 2 -y 2 with q = (7r,0,7r) for AE = —0.5 eV, in 
the BaNi2As2 system. Since AE = corresponds to the 
tetragonal phase, the AF Oxz(yz) fluctuations would be 
dominant in the tetragonal phase. Thus, the supercon- 
ductivity in the tetragonal phase with T c ~ 3 K recently 
reported by Kudo et al£& might be caused by the AF 
Oxz(yz) fluctuations with q = (7r,7r,7r). In contrast, a 
large energy level shift AE would be required to real- 
ize the AF Ox 2 -y 2 fluctuations in the tetragonal phase 
(The necessity of the energy level shift AE means that 
the structure transition is of first order). If we could re- 
alize the second-order structure transition in some way 
such as by applying uniaxial pressure, the AF Ox 2 -y 2 
fluctuations may be stronger, and higher- temperature su- 
perconductivity might be realized in Ni-based supercon- 
ductors. 



VI. SUMMARY AND FUTURE PROBLEMS 

To summarize, we investigated the electronic state and 
structure transition of BaNi2As2- The Coulomb and 
quadrupole-quadrupole interactions were treated within 
the RPA method. The former interaction is not impor- 
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tant since U c is very large and, as a result, spin fluctua- 
tions are very small. Owing to the latter interaction, two 
types of quadrupole fluctuations develop in BaNi2As2: 
AF Oxziyz) an d AF Ox 2 -y 2 fluctuations. The lat- 
ter fluctuations are the origin of the experimental zigzag 
chain formation, 27 and the former fluctuations would 
then be the origin of the strong-coupling superconductiv- 
ity observed in BaNi2(Asi_ x P x ) 2 £&T c m Ni-based super- 
conductors would increase further by developing stronger 
AFQ fluctuations. We conclude that both Ni- and Fe- 
based superconductors are characterized as strong orbital 
fluctuating metals, although origins of their orbital fluc- 
tuations are different. 

In the present paper, we discussed the zigzag chain for- 
mation due to the AF Ox 2 -y 2 fluctuations for AE = 
—0.5 eV, and the resultant crystal structure is or- 
thorhombic. Experimentally, both longer bond (~ 3.1 A) 
and shorter bond (~ 2.8 A) are slightly modulated in a 
staggered way. 27 and the crystal structure becomes tri- 
clinia The origin of this triclinic structure formation is 
an important future problem. 

We also comment on the role of VC: In Ref.[i|, we found 
that the Aslamazov-Larkin-type VC due to Coulomb in- 
teraction causes strong orbital fluctuations, which would 
be the main origin of the orthorhombic phase transi- 



tion as well as superconductivity in Fe-based supercon- 
ductors. However, this mechanism will be unimportant 
in BaNi2As2 because of the large U c . On the other 
hand, the VC due to electron-phonon interaction might 
be important in BaNi2As2- By this mechanism, the 
AF Or fluctuations give rise to the ferro Ov order if 
Tr{O r Or'} ^ (two-orbiton term)££ In the present 
model, this mechanism can realize the ferro Osz 2 -r 2 or- 
der for AE = —0.5 eV, which causes the change in the 
inter-layer distance with preserving the tetragonal sym- 
metry. It might be related to the abrupt change in c-axis 
resistivity at T s . 27 This is another important future prob- 
lem. 
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